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Abstract 

We investigate the extension of seh-injecting laser wakeheld experiments to the regime that wih 
be accessible with the next generation of petawatt class ultra-short pulse laser systems. Using 
linear scalings, current experimental trends and numerical simulations we determine the optimal 
laser and target parameters, i.e. focusing geometry, plasma density and target length, that are 
required to increase the electron beam energy (to > 1 GeV) without the use of external guiding 
structures. 
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I. INTRODUCTION 



The field of laser driven electron acceleration has attracted significant international at- 
tention since the first observations of quasi-monoenergetic electron beams [H 121 E] ■ Results 
on the production of mono-energetic beams in self-injected laser wakefields from both exper- 
iments (solid symbols) and particle-in-cell simulations [ll[5l|6l[71[8l|9l[ini[IIl[12l[13l[Tl[T5] 
have been collated in figure [T] 

The maximum electron energy and acceleration length required in a laser wakefield ac- 
celerator can be derived from simple considerations. The maximum length over which a 
wakefield can accelerate particles is typically governed by the dephasing length, that is the 
length over which a relativistic {v — > c) electron will overtake the plasma wave, which travels 
at the group velocity of the laser pulse {vg ~ c(l — rie/ric)^^^), by half the wavelength of 
a relativistic plasma wave (Ap = 2Tic/ujp). Ue is the plasma density, ric = meoUQ/e"^ is the 
critical plasma density for electromagnetic wave propagation, where Uq is the laser frequency 
{ric = 1.75 X 10^^ cm~^ for 800 nm radiation). The dephasing length for rif. ^ is given 
by the formula 

The maximum electron energy gain that a wakefield accelerator can produce can be 
estimated by integrating the electric field over a dephasing length (assuming a sinusoidal 
and non-evolving electric field). The electric field amplitude required to trap electrons that 
are initially at rest is given by Eq = mcUp/e. Note that this is not the field at which ID 
wave-breaking occurs, -Ewb = (2(7,/, — l))^/^-E'o[16j, where 7,^ is the Lorentz factor associated 
with the phase velocity of the plasma wave {v^ ~ Vg). -Ewb is the electric field at which 
the plasma wave can trap and accelerate an electron initially moving at —v^ to one moving 
at +v^. The field required to trap an electron at rest is more relevant to the transverse 
injection that is important for narrow energy spread beams. This is because the electrons 
that are travel across the back of the first wave period, or 'bubble', have low longitudinal 
momentum and can therefore be trapped at lower field strengths, typically on the order of 
Eq, than the ID limit. 

This leads to an expression for the maximum energy an electron can gain from a plasma 
wave of amplitude Eq. 

Wmax = 2mc'— (2) 
2 




FIG. 1: Collation of reported data from various experiments and 3D particle- in-cell simulations. 
Closed circles: electron beam energy from experiments. Open circles: electron beam energy from 
simulations. Closed squares: interaction length from experiments. Open squares: interaction 
length from simulations. The solid line corresponds to equation [2] and the dashed line corresponds 
to equation [1} Data taken from [IlElElllElElITlElElIinKIIKniIISlEKIS] 



Equations [T] and [2] are also plotted on figure 1. It should also be noted that this figure 
is on a log-log scale, so that any scatter from the theoretical curve appears diminished. 
Nevertheless, the observed electron energies and interaction lengths appear to follow these 
simple scaling laws over a range of plasma densities. It should be noted that not all the 
data points are for similar laser systems. Data has been included from experiments with 
laser power as low as 2 TW and as high as 40 TW. Published 3D Simulation data has been 
included for laser powers over the range 10 - 300 TW. Both experimental and simulation 
data presented includes self-guiding and external guiding channel results (from less than 1 
mm to > 10 mm). 

This collection of data clearly shows the general trend that to increase the electron beam 
energy experiments must move to lower plasma density and longer interaction lengths. Stable 
0.5 GeV acceleration has been achieved at LBNL following this approach, with some shots 
reaching the GeV level [9] . It should be noted that more detailed scalings were proposed in 
[TTj and [H] that include the dependence on the laser intensity but also assume the pulse 
dimensions are matched to the bubble or plasma wave dimensions which is not always the 
case in the experimental data presented. These experiments typically rely on significant 
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pulse evolution before electron injection can occur. 

It is not sufficient simply to reduce the plasma density to increase the electron beam en- 
ergy; the laser parameters must also be modified so as to reach this mono-energetic regime. 
The self-similar behaviour of these accelerators with density requires that the pulse dimen- 
sions scale with \fn^ while maintaining a minimum intensity; this requires the laser energy 
to increase with decreasing density. Tsung et al [3] and Lu et al [TU [19] have shown from 
simulations and analytical theory, that for injection at the back of the first wave period, a 
minimum intensity threshold of approximately Oq > 3 is required. 

While pulse evolution (i.e. self-focusing and pulse compression) has hitherto played a 
crucial role in reaching self-injection with 10-100 TW lasers, it has also been attributed to 
some of the remaining shot-to-shot variability of the electron beam parameters. Experimen- 
tal studies have shown that the stability of the electron beam is increased when the pulse 
and plasma parameters are chosen such that the beam waist Wq is matched to Ap, minimis- 
ing self-focusing effects [20] . It has also been experimentally verified that, in this regime, 
self-focusing tends to produce exit mode profiles with a beam waist approximately equal to 
the plasma wavelength [21j for laser powers above the critical power for self-focusing. 2D 
simulations also show that, for a range of plasma densities and focal geometries, self- focusing 
reduces the beam size until it is approximately wq ~ Ap after which stable propagation oc- 
curs. Fig. |2] shows experimental and 2D Osiris simulation results that have been found to 
support this statement with 15 TW laser systems. 

A number of laser systems are now being constructed with laser powers 100 TW - IPW 
and with laser pulse durations < 100 fs. One of the main goals of such lasers is their use as 
drivers for next generation particle acceleration experiments. For example, a single beam of 
the Astra Gemini laser, which is under construction at the Rutherford Appleton Laboratory, 
is expected to produce a laser pulse of duration ^ 30 fs with pulse energy of £^ ~ 15 J on 
target [22] ■ 

To calculate the minimum density at which a laser such as Astra Gemini will be able to 
self-trap electrons we use a simple model [20]. We assume that the beam self-focuses to a 
matched spot size, wq ~ Ap and require that the minimum intensity of the evolved pulse 
required to self-trap electrons corresponds to a normalised vector potential of ^ 3.2. This 
value of at was taken from experimental observations of the trapping threshold on 20 TW 
experiments [2D] and in agreement with the simulations reported in [5]. The relationship 
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FIG. 2: Results from a recent Astra 15 TW experiment (black squares) and 2D PIC simulations 
(grey circles) showing the dependence on the laser spot size (FWHM intensity « 1/e^ radius) with 
plasma density (experimental data taken from [21]). The solid line indicates Xp = 27rc/wp. The 
dashed line shows the density below which the laser in the experiment was below the critical power 
for self-focusing Pc/[GW] = 17.3(nc/ne). 



between the laser pulse energy £l, pulse duration (FWHM) tl, and plasma density Ue can 
be found using by comparing the power of the laser pulse: 

~ ° ° ^ 

and the critical power for self-focusing: 



which produces 



Pc 8 A2 



Inserting the threshold vector potential required for trapping, at and using the fact that 
the peak laser power Pl 0.9Sl/tl for a gaussian temporal profile, and that the spot size 
evolves towards wq ~ Ap yields the result 



To verify this simple model we calculated the minimum laser energy required to achieve 
trapping for a range of reported experiments and simulations [Il[2l[3llll|5l[6l[7l[8l|9|, [TOj [TTl 
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FIG. 3: Calculated minimum laser pulse energy required for self-trapping from equation |6] for 
reported experimental and simulation data against the actual laser energy. Solid circles: self- 
guided experiments. Solid squares: self-guided simulations. Open circles: guided experiments. 
Open squares: guided simulations. The dashed line represents £l = £t calculated for at = 3.2. 
Data taken from [Ilia|3lia|5llS[3[8lllll0llIIl[l2lll3l[Il[l5l|20l|2a 



[I3[T3l[Tl[T5l|20l|23l[2l|25l[26l|271|2^ The calculated energy threshold expression was 
modified to include the effect of pulse compression by replacing in equation[6]by /ujp (i.e. 
ctl = Ap/2). In one case we used the measured pulse duration after compression [21]. In 
figure [3] we plot the calculated minimum energy required for self-trapping against the actual 
energy in the laser pulse. Note that only one self-guided experimental result [10| and no 
self-guided simulation results have a laser energy lower than the predicted threshold. This 
implies that the energy threshold model is reasonably accurate for self-guided experiments. 
Apart from one experimental self-guided data point, all the other points that have Sl < St 
are experiments or simulations where pre-formed guiding structures were used. This could 
indicate that the trapping mechanism in some guided experiments is not simply self-trapping 
or it could indicate that the presence of a radial density profile lowers the required intensity 
in the self-focused and compressed pulse as discussed in [T8] . 

The minimum density at which we expect trapping to occur (from equation [6] ) is shown 
as a function of laser energy for = 30 fs in figure |4j For Sl = W J we therefore expect 
the minimum density at which self-injection will be possible will be ne ~ 1 x 10^^ cm~^ and 
that the corresponding electron energy will be W^ax ~ 1.5 — 2 GeV. The dephasing length 
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FIG. 4: Minimum density at which trapping wih occur as a function of laser energy (left axis) and 
the expected eletron beam energy (right axis). The threshold was calculated using equation [6] with 
r = 30 fs and at = 3.2. The dashed line indicates the position of a iSl = 10 J laser (e.g. Astra 
Gemini) . 

at this density from equation [l] is close to 6 cm, which would require significant self-guiding. 

II. NUMERICAL MODELLING OF GEMINI ELECTRON ACCELERATION EX- 
PERIMENTS 

A series of simulations investigating potential electron acceleration experiments have been 
carried out for a = 30 fs, = 10 J A^, = 800 nm laser system. The simulations were 
performed using the particle in cell code Osiris in 2D3V slab geometry. In these runs the 
laser propagates in the x direction, the slab lies in the x — y plane and the fields and particle 
momenta have components in x, y and z. The laser was polarized in the x — y plane. We use 
stationary ions. The simulations were performed on the 48 node 'Caesar' cluster and on the 
CXI supercomputer at Imperial College. The simulation resolution was carefully chosen to 
minimise numerical dispersion errors while maintaining an acceptable run size. Typical run 
parameters were Aa; = 0.2 c/ujq ; Ay = 0.8 c/t^o ; At = 0.199 l/tuo in a simulation box size 
of up to (1600 X 1600) c/ujq ~ 200 x 200 /im. The simulation box moves in the direction of 
the laser propagation at the speed of light. Runs were performed for propagation distances 
as large as 1 cm. The physical parameters of the runs presented in this report are shown in 
table E 
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Simulation density Beam waist Rayleigh Range vector 



1 

number 


/ -3 
/ cm 


Wo / /im 


Zr I mm 


4- 4- ■ 1 

potential oq 


1 


1.05 X 10^® 


5 


0.1 


18.8 


2 


1.05 X lO^s 


10 


0.4 


9.6 


3 


1.05 X 10^^ 


20 


1.6 


4.8 


4 


1.05 X 10^^ 


34 


4.5 


2.8 


5 


2.10 X lO^s 


20 


1.6 


4.8 



TABLE I: Physical parameters used for the simulations presented, ti = 30 fs (FWHM) and 
El = ^ for all the simulations. 

The first 4 simulations were performed at tlie minimum plasma density for self-trapping, 
chosen using equation |6} of rig = 1 x 10^^ cm^'^ but varying the spot size. In experiments 
the focusing geometry is not a simple parameter to change, and in PW class systems the 
necessary increase in beam size (~ 10 cm diameter) to prevent damage to optics increases 
the size and cost of vacuum systems. Long focal lengths further exacerbate this situation 
and it is therefore advantageous to study the minimum focal spot size at which self-guiding 
will occur at the expected minimum density for self-trapping. 

Figure [5] shows the evolution of the laser electric field envelope in simulations 1-4. Each 
vertical slice in the image corresponds to the transverse profile of the laser envelope after 
integration along the propagation direction, x. For small spot sizes, Wq = 5 and 10 /im, 
which correspond to wq < Ap, the pulse clearly diverges and is not significantly self-guided. 
Indeed these simulations were halted after 2 and 4 mm propagation due to interaction of 
the diffracting pulse with the box boundary. As the spot size becomes close to the plasma 
wavelength guiding is observed. For wq = 20 and 34 /xm there is still significant laser 
intensity (ao ~ 2 — 3) after a propagation distance of 1 cm. The Wq = 34 /im simulation has 
guided over approximately 2Zji {Zr = ttWq/Al is the Rayleigh range) before the simulation 
was stopped. However the wq = 20 fim case has guided for over 6Zr. 

Self-guiding alone is insufficient to produce an electron beam. The self-guided pulse 
must maintain a sufficiently high intensity to produce self-injection. In the following we 
concentrate on self- injection in the first wave period since 2D simulations can over estimate 
self-injection in the trailing periods compared with full 3D runs [31J. In fig. [6] we show the 
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FIG. 5: Evolution of the transverse laser pulse envelope at a plasma density of ng = 1 x 10 cm~ 
for various beam waists. Simulations 1-4 are shown. The transverse size of simulations 3 and 4 
was increased to 200 /xm compared to 100 fim for simulations 1 and 2. 

electron density distribution of the first plasma wave period from three simulations (numbers 
4,3 and 5 in table 1) after the laser has propagated for 9.6 mm. In these plots the laser 
propagation direction is left to right. 

Simulation 4, at a density of ne = 1 x 10^^ cm~^ and a spot size matched to the plasma 
wavelength Ap has not injected in the first period even at this late stage. In simulation 4, 
ao = 2.8, which is close to, but just below, the expected threshold for injection {at ~ 3). 
Pulse modification would be required before injection could occur, however as the pulse 
waist is already matched to Ap and the pulse length (cr = 9 /xm) is less than Ap/2 we might 
expect minimal self-focusing and compression. 

Simulation 3, also at a density of ne = 1 x 10^^ cm~^ has also failed to produce self- 
injection in the first period. In this case oq = 4.8 which is above the expected threshold. 
This failure is explained by the fact that in simulation 3 wq = 20 fim, which is less than 
Ap = 34 /im so that some diffraction occurs initially, reducing the intensity before self-guiding 
occurs and so the plasma wave does not reach a sufficiently large amplitude. 

Simulation 5 shows the effect of moving to a higher plasma density {rie = 2 x 10^* 
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FIG. 6: Electron density profile after 9.6 mm propagation from three simulations investigating 
injection threshold. The laser travels from left to right. 

cm~^), while maintaining the spot size of simulation 3. In this case the plasma wavelength 
is Xp = 23 yum, which is very close to the spot size wg = 20 /im, and hence sufficient 
laser energy is trapped in the first plasma wave period and the intensity was high enough 
to cause self-injection after 4.2 mm propagation (i.e. some pulse evolution was required 
before injection occurred). The electron density plot shown in figure |6] is after 9.6 mm of 
propagation, and therefore shows an electron bunch that has been accelerated for over 5 mm. 
The fact that some pulse evolution was required prior to injection may indicate that the 
most stable beams from such a laser system may be achieved at a slightly higher plasma 
density, and therefore with a slightly lower electron beam energy. 

Fig. [7] shows the accelerating field at the point where injection has just occurred. The 
peak accelerating field is close to = 5 mcup/e = 0.7 GeV/mm. We therefore expect that, 
if this field were maintained over the entire acceleration length of 5 mm then electron energies 
close to 3 GeV would be obtained. This electric field at the back of the bubble is above 
the trapping threshold for stationary electrons, Eq, but is significantly less than the ID cold 
wavebreaking field for this plasma density, -Ewb ~ 40i?o, indicating that multi-dimensional 
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FIG. 7: Longitudinal electric field {E^) of the first plasma wave period in simulation 5 after 4.2 
mm propagation. Blue (red) regions correspond to accelerating (decelerating) fields. 



effects are significant in determining the exact electron trajectories. That the field in the 
bubble is significantly larger than Eq is somewhat at odds with the scaling presented in 
fig. [Tj Using that scaling we would predict peak electron energies of 0.85 GeV at the plasma 
density of simulation 5 {rie = 2 x 10^* cm~^). 

An electron spectrum obtained from simulation 5, considering only electrons that would 
pass through an electron spectrometer with a 25 mrad acceptance cone (corresponding ap- 
proximately to a typical electron spectrometer) is shown in fig. [s] This shows a quasi- 
mononenergetic spectrum at approximately 2 GeV (2 % relative energy spread) indicating 
that the size accelerating field does decrease slightly over time, possibly due to laser pump 
depletion or beam loading effects. We note that in this case the maximum energy is close 
to Wmax ~ 2ao'^c^('^c/'^e), which is consistent with 3D non-linear scalings including the 
dependence on oq [IS]- 



III. SUMMARY 



In a series of particle-in-cell simulations we have investigated the parameters of interest 
for the next generation of experiments on self-guided laser wakefield acceleration. Using 
simple scaling laws we identified the regions of interest and have then verified and refined 
these using 2D PIC simulations. We have considered the appropriate focusing geometry and 
plasma target requirements necessary to produce quasi-monoenergetic electron beams to the 
multi-GeV level. The simulations show that a focal spot size on the order of 20 /im (i.e. 
//20 focusing) is ideal and that plasma targets of centimetre length capable of producing 
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FIG. 8: Electron energy spectrum (number of electrons per relative energy spread) after 9.6 mm 
propagation in simulation 5. The electron spectrum is calculated for electrons travelling within a 
cone angle of 25 mrad. The FWHM energy width of this beam is 2%. 

densities up to a few 10^^ cm~^ are suitable. This indicates that multi-GeV beams from a 
self-guided wakefield accelerator should be experimentally realisable in the near future. 
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